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ABSTRACT 

The configuration and evolution of the magnetic field in star-forming cores are investigated in order 
to directly compare simulations and observations. Wc prepare four different initial clouds having dif- 
ferent magnetic field strengths and rotation rates, in which magnetic field lines are aligned/misaligned 
with the rotation axis. First, we calculate the evolution of such clouds from the prestellar stage until 
long after protostar formation. Then, we calculate the polarization of thermal dust emission expected 
from the simulation data. We create polarization maps with arbitrary viewing angles and compare 
them with observations. Using this procedure, we confirmed that the polarization distribution pro- 
jected on the celestial plane strongly depends on the viewing angle of the cloud. Thus, by comparing 
the observations with the polarization map predicted by the simulations, we can roughly determine 
the angle between the direction of the global magnetic field and the line of sight. The configuration of 
the polarization vectors also depends on the viewing angle. We find that an hourglass configuration 
of magnetic field lines is not always realized in a collapsing cloud when the global magnetic field is 
misaligned with the cloud rotation axis. Depending on the viewing angle, an 5"-shaped configuration of 
the magnetic field (or the polarization vectors) appears early in the protostcllar accretion phase. This 
indicates that not only the magnetic field but also the cloud rotation affects the dynamical evolution 
of such a cloud. In addition, by comparing the simulated polarization with actual observations, wc 
can estimate properties of the host cloud such as the evolutionary stage, magnetic field strength, and 
rotation rate. 

Subject headings: ISM: clouds — ISM: magnetic fields — methods: numerical — polarization — stars: 
formation 



1. INTRODUCTION 

The magnetic field is believed to play a key role in the 
star formation process. Stars form in molecular cloud 
cores. In the collapsing cloud core, the magnetic field 
transfers excess angular momentum and promotes fur- 
ther gas contraction to form a protostar. After the pro- 
tostar formation epoch, the protostellar outflow appears 
to be driven by the magnetic force. The protostcllar 
outflow influences the interstellar medium and seems to 
promote subsequent star formation. 

To investigate the effects of the magnetic field, wc 
need to determine its strength and morphology. The 
strength of the magnetic field in molecular clouds can 
be determined by Zeeman splitting observations. Many 
observations have shown that in a molecular cloud core, 
the magn etic energy is a s large as the gravitational en- 
ergy (e.g.. ICrutcherlll999l ). Thus, the magnetic field in a 
gravitationally contracting cloud cannot be ignored. On 
the other hand, the morphology of the magnetic field in 
molecular clouds has been obtained by polarization mea- 
surements of thermal dust emission. 5 In past studies, 
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J Because the dust major axis is aligned perpendicular to the in- 
terstellar magnetic field, thermal emission from such dust is linearly 
polarized, and the B-vector of the electromagnetic wave indicates 



the direction of the magnetic field lines at the molecu- 
lar cloud scale was estimated from the polarization vec- 
tors (e. g..lMoneti et al.lll984lTamura fc Satol ll989D. Re- 
cently, iGirart et al.l (|2006l I2009TT obtained a fine struc- 
ture of collapsing cloud cores by high-angular-resolution 
measurements of polarized emission. They showed a very 
clear hourglass structure of magnetic field lines (or po- 
larization vectors). The hourglass structure is believed 
to be an evidence of gravitational contraction, because 
the magnetic field, which is well coupled with neutral 
gas, is pinched toward the center of the gravitationally 
collapsing gas cloud. 

To understand the effects of the magnetic field on the 
star formation process, information about both the mag- 
netic field strength and its morphology are desirable. 
However, it is difficult to determine the magnetic field 
strength at a small scale (or in a high-density gas re- 
gion) by resolving the internal structure of the collapsing 
cloud, which can only be done thro ugh Zeeman splitting 
observations (jCrutcher et al.|[2~010[ ). 

On the other hand, observation by dust polariza- 
tion measurement can resolve small-scale structure with 
(near-future) interferometers such as the Submillime- 
tcr Array (SMA) and the Atacama Large Millimc- 
ter/submillimctcr Array (ALMA). Thus, we can extract 
various types of information on the star formation pro- 
cess from polarization observations. 

Polarization observations provide information of the 
magnetic field morphology projected on the celestial 
plane. Because, in reality, magnetic field lines have a 
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three-dimensional configuration, we need to reconstruct 
a three-dimensional configuration of magnetic field lines 
from the observed polarization vectors, considering the 
projection effect. However, such approach is very com- 
plicated, and some parameters, such as the projection 
angle, are uncertain. Instead, in previous studies, re- 
searchers develop ed models of the magnetic field lines 
witho ut rotation ([Goncalves et al.l 120 08; Pad ovani et al.l 
120121) or with the effect of rotation fe.g. lFrau et al.ll2011l ). 
Then, they compared their model with observations, in 
which the projection effect of the magnetic field lines are 
taken into account. However, since magnetic field lines 
modeled in such studies are somewhat idealized, it is dif- 
ficult to investigate effects of different cloud conditions 
and time evolution on magnetic field lines. Thus, al- 
though such models may roughly explain the large-scale 
(or molecular cloud scale) magnetic field structure, we 
need a more sophisticated model based on the cloud 
collapse simulation to explain recent (or future) high- 
resolution observations. 

Recent observations could resolve the polarization dis- 
tribution even around the region close to the protostar, 
where the effect of cloud (or disk) rotation, protostellar 
outflow, and the circumstellar disk is expected to be sig- 
nificant, in order to investigate (or model) the magnetic 
field lines. Thus, we need to construct a more realistic 
magnetic field line model considering these effects. 

The evolution of magnetized clouds has recently 
been investigated by magnctohydrodynamic (MHD) 
simulations. Suc h simulat i ons reproduced the pro- 
toste l lar outflow (ITomisakal 120021: iBaneriee fc Pudritzl 



2006t iMachida et alJl2007HHennebelle fe Fromaiidl2008t 



Puffin fe Pudritzll2009D ~ and demonstrated the formation 



of the circumstellar disk in the c ollapsing magnetized 
cloud dHennebelle fc Ciardil 120091: IMachida et all 120111 : 
iDuffin et al.l 1201 lh . These studies also showed the mag- 
netic field lines in three dimensions or the magnetic field 
vector on an arbitrary cutting plane. The magnetic field 
obtained from the numerical simulations is considered to 
be more realistic than a simple theoretical model. How- 
ever, we cannot directly compare the magnetic field ob- 
tained from numerical simulations with polarization ob- 
servations without introducing the effect of projection 
and thermal dust emission. Thus, we need a further at- 
tempt to compare numerical simulations with observa- 
tion results. 

To compare n umerical simulations with observations, 
ITomisakal (|2011| ) calculated the polarization of thermal 
dust radiation from two-dimensional axisymmctric MHD 
simulation results, assuming a constant emissivity per 
unit mass. He called this procedure "observational vi- 
sualization." In this study, first, we calculate the cloud 
evolution from the prestellar cloud core stage until long 
after protostar form ation. Th e n, acc ording to the pro- 
cedure described in ITomisakal (|2011l) , we calculate the 
polarization distribution predicted by MHD simulations 
and discuss the relationship between the polarization 
vectors and the evolution ary stage of the cloud or pro- 
tostar. ITomisakal ()2011l ) assumed magnetic field lines 
parallel to the rotation axis in the initial cloud because 
he used the data obtained from axisymmctric simula- 
tions. In our study, we investigate both aligned and mis- 
aligned magnetic fields. The initial magnetic field lines 
are aligned with the initial rotation axis in the former 



case, whereas they are inclined from the rotation axis 
in the lat ter case. Part of our results was already pre- 
sented in iShinnaga et al.l (|2012f) , in which we compared 
the polarization vectors derived from numerical simula- 
tions with polarization observations in a massive star- 
forming region and determined various cloud parameters 
and the evolutionary stage of the protostar. This pa- 
per is structured as follows. The framework of the MHD 
simulation is given in §2. The numerical method for cal- 
culating the polarization is presented in §3. The results 
of the MHD simulation and polarization calculation are 
shown in §4. We discuss the configuration of magnetic 
field lines and compare our results with observations in 
§5. We summarize our results in §6. 

2. MODEL SETTING FOR CLOUD COLLAPSE 
SIMULATION 

2.1. Basic Equations and Initial Settings 

To study the evolution of a magnetized cloud core, we 
solve the three-dimensional resistive MHD equations, in- 
cluding self-gravity: 



(i) 



v + p ( v . v)« = -VP - — B x (V x B) - P V#2) 
at 4-7T 



— = V x (v x B) + ??V 2 B, 

V 2 (f) = 4nGp, 



(3) 
(4) 

where p, v, P, B, rj, and <f> denote the density, velocity, 
pressure, magnetic flux density, resistivity, and gravita- 
tional potential, respectively As the gas pressure, we 
adopt the isothermal equation of state as 



P = C s n P, 



(5) 



where c s .o = 190 ms -1 . The isothermal approx- 
imation is justified i n the low-density gas region. 
iMasunaga fc Inutsukal (|2000f ) showed that in the col- 
lapsing cloud, the gas behaves isothermally for n <J 

10 10 cm -3 and adiabatically for n ;> 10 10 cm~ 3 . In 
this study, we calculated only a low-density gas region 
(n < 10 10 cm~ 3 ) using a sink cell treatment (see W2.2\i . 
To safely calculate the gas evolution around or inside the 
sink, we adopted an artificially large resistivity in equa- 
tion ([3]), as described in §2.21 

As the initial state, we assume a spherical cloud with a 
critical Bonnor-Ebcrt (BE) density profile Pbe, in which 
a uniform density is adopted outside the sphere (r > R c , 
where R c is the critical BE radius). For the BE density 
profile, we adopt an isothermal temperature of T = 10 
K and a central number density of hq = 5 x 10 5 cm -3 . 
For these parameters, the critical BE sphere has a radius 
of R c = 6.5 x 10 3 AU and a mass of 1.1 M . The gravi- 
tational force is ignored outside the host cloud (r > R c ) 
to mimic a stationary interstellar medium. Thus, only 
the gas inside r < R c can collapse to form the protostar. 
Because the critical BE sphere is in an equilibrium state, 
we increase the density by a factor of / = 1.67 to pro- 
mote contraction, where / is the density enhancement 
factor. 

In the initial cloud, the uniform magnetic field Bq par- 
allel to the z-axis is adopted in the entire computational 
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TABLE 1 
Model parameters 
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domain. The cloud rotates rigidly and the rotation axis 
is inclined with respect to the initial magnetic field lines. 
Thus, the initial angular velocity is given by 

/ sin 8 \ 

= n 1 o , (6) 

\COS SqJ 

where S is the angle between the initial magnetic field 
lines (or z-axis) and the rotation axis. 

The initial cloud is characterized by three parameters, 
Bq, S7o, and So, which are summarized in the first to third 
columns of Table 12.11 The non-dimensional quantities 
ao , /3o , and 70 are also listed in the forth to sixth columns 
of Table 12. 1[ where ao, (3o, and 70 are the ratios of the 
thermal, rotational, and magnetic energies, respectively, 
to the gravitational energy inside the initial cloud (R < 
Rc)- 

We estimate the mass-to-flux ratio of the initial cloud 



M 



M 



$ nR 2 cl Bo ' 



(7) 



where M is the mass contained within the critical ra- 
dius R c and $ is the magnetic flux threading the initial 
cloud. A critical value of M/$ exists below which a cloud 
is supported against self-gravity by the magnetic field. 
For a cloud with uniform density, iMouschovias fc Spitzerl 
(|1976t) derived a critical mass-to-flux ratio, 



$ A 



3tt 



1/2 



(8) 



where the constant ( = 0.53 for unifor m spheres (C = 
0.48 a ccording to careful calculation by iTomisaka et al.l 
[T988aHbT ). For convenience, we use the mass-to-flux ratio 
normalized by the critical value as 



AT 



M 



(9) 



The models have the normalized mass-to-flux ratio /i in 
the range of 2.3 < /x < 12.6. Thus, the initial clouds are 
magnetically supercritical. T his is supported by recent 
observations (jCrutcherl 119991 ) . The normalized mass-to- 
flux ratio for each model is also listed in the seventh 
column of Table HHJ 

2.2. Numerical Method and Sink Cells 

To calculate on a large spatial scale, t he nested grid 
metho d is adopted (for details, sec Machid a et al.l 
l2005allbT ). Each level of a rectangular grid has the same 
number of cells (64 x 64 x 64). The calculation is first 
performed with five grid levels (1 = 1 — 5). The box size 



of the coarsest grid, I = 1, is set to 2 5 i? c . A new finer 
;rid is generated bef ore the Jeans condition is violated 
Truelove et al|[l997l ). The box size of the first level of 
grid (I = 1) is L\ ~ 2.1 x 10 5 AU, whereas the maximum 
grid level, Z max = 9, has a box size of Lg = 820 AU and 
a cell width of A i=9 = 13 AU. 

With these settings, we calculated the cloud evolu- 
tion until a large fraction of the total cloud mass fell 
onto the protostar. To perform long-term calculation of 
the collapsing cloud, we adopted a sink at the center of 
the cloud. We started the calculation without a sink 
and calculated the cloud evolution for the prestellar gas- 
collapsing phase without the sink. Later, we identified 
protostar formation in the collapsing cloud core when 
the number density at the cloud center exceeds n > n t hr , 
where n t hr is the threshold density for our model. After 
protostar formation, wc calculated the cloud evolution 
with the sink. 

To model the protostar, we adopted a fixed sink with 
a radius of r s ; n k = 15 AU composed of sink cells only 
around the center of the computational domain. Because 
we added no non-axisymmetric perturbation to the initial 
state, the protostar (or the center of gravity) does not 
move and remains at the center of the computational 
domain during the calculation. In the region r < r s i n ]< = 
15 AU, gas having a number density of n > nthr = 3 x 
10 9 cm -3 (pthr = 1-2 x 10~ 14 gcm~ 3 ) is removed from 
the computational domain and added to the protostar 
as a gravitating mass in e ach timestep (for details, see 
IMachida et~aLll2009l [2010h . Thus, for each timestep, the 
mass accreted onto the protostar is calculated as 



AM, 



\p{h3, k ) - (hbx) dV. 



(10) 



r<f Bink 



To avoid artificial amplification of the magnetic field 
around or inside the sink, we adopted the resistivity in 
equation ([3]) as 



V 



740 / T 



10K 



2 — i 
em's , 



(11) 



where T (10 K) is the gas temperature and X e is the 
ionization degree of the gas, which is expressed as 



X. 



5.7 x 10" 



lem 



(12) 



Equation (TTT1) coincides with the resistivity of th e col- 
lapsing gas cloud de rived in iNakano et al.l (|2002l ) and 
IMachida et al.l ()2007l ) with a numerical factor c v = 1, in 
which the magnetic Reynolds number Re becomes Re < 

I at n ~ 10 12 cm -3 , and the magnetic field begins to 
dissipate by Ohmic dissipation for n > 10 10 — 10 12 cm" 3 
(for details, see IMachida et al.ll2007l l2008al |bl). 

The magnetic field is highly amplified around the pro- 
tostar without the magnetic dissipation because the mag- 
netic flux continues to flow into the sink cell. Such 
field may induce artificial effects in the region outside 
the sink cell. To avoid this, we chose c v = 10 to dissi- 
pate the magnetic field inside or around the sink; conse- 
quently, the magnetic Reynolds number becomes Re < 1 
at ?i ~ 10 11 cm -3 and the magnetic field dissipates for 

II ;> 10 9 cm~ 3 . In our setting, the gas continues to 
flow into the sink cell, while the magnetic field dissi- 
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pates around and inside the sink. We confirmed that 
this treatment (i.e., artificial dissipation of the magnetic 
field around the sink) has little effect on the outer gas- 
collapsing region (see also £I5.4|) . Therefore, we can safely 
calculate the cloud evolution for a long period. 

3. NUMERICAL METHOD FOR POLARIZATION 
CALCULATION 

After the MHD simulations were conducted,we calcu- 
lated the polarization of thermal dust emission using the 
results of MHD simul ations according to the formulae in 
iTomisakal (pOll (see iFiege fe Pudritzl [2000h . Here, we 
briefly summarize the formulation. 

The purpose of this study is to calculate the polar- 
ization vectors and polarization degrees of thermal dust 
emission. For simplicity, we assume that the entire re- 
gion in the cloud is optically thin and isothermal. We 
define the observed direction as the direction of a unit 
vector n. The observational grid normal to n has the 
following horizontal (£) and vertical (77) axes: 



i e ^ e '- n |" , (e z -n^0) 

e z — (e z -n)n ' V z ' ' 

-e x , (e z • n = 0) 



(13) 



and 



= e,, x n. (14) 
Then, we calculated the relative Stokes parameters, 



q = I p cos 2ip cos 7 ds 



and 



u = I p sin 2tp cos 2 7 ds, 



(15) 



(16) 



where the integration is performed along the line of sight 
n, and the angles ip and 7 are between the projected 
magnetic field direction and the rj-axis and between the 
magnetic field direction and the £77 plane, respectively. 
The polarization degree is given as 



P 



Pa- 



S - 7J0S2 



(17) 



from two quantities obtained by integration of the den- 
sity along the line of sight as 



and 



£2 = 



£ = / pds 



P 



cos 7 1 . 
— — ds. 



(18) 



(19) 



In equation (fTT)) . the numerical factor po is set to 0.15 
in order to fit the observational maximum polarization 
degree in interstellar space. 

Using equations (fT5"]) -(fT9 j) . we calculated the polariza- 
tion at various viewing angles for different datascts taken 
from simulations. The directions of the viewing angles 
are specified by the angles 9 and 4>: 9 represents the in- 
clination angle and <j> is the azimuthal angle. The case of 
9 = 0° corresponds to the pole-on view, whereas 9 = 90° 



simulation grid (x,y,z) 

celestial plane (f , T\ ) 




observer 



Fig. 1. — Relationship between simulation grid and observation 
grid. 

corresponds to the edge-on view. Because of symme- 
try, we examine the inclination angle in the range of 
0° < 9 < 90°. The relationship between the simula- 
tion grid (xyz coordinate system) and the observation 
grid (£17 coordinate system) is shown in Figure [TJ 

4. RESULTS 

In this section, we show the results of MHD simulations 
in §4.11 and their observational visualizations in §4.21 for 
the four models (models 1-4) listed in Table 12.11 We use 
different coordinate systems in N4. II and £14.21 because we 
observe the collapsing cloud from an arbitrary direction 
in §L2l 

The configurations and directions of the magnetic field 
lines and rotation axes for the models are schematically 
illustrated in Figure [2 For all the models, we adopted a 
uniform magnetic field that has the same direction but 
different strengths for each model. As seen in the fig- 
ure, model 1, which is the standard model, has a uni- 
form magnetic field parallel to the z-axis without rota- 
tion. Model 2 has a magnetic field of the same strength 
and configuration but with rotation; the rotation axis is 
parallel to the magnetic field lines (or parallel to the z- 
axis). On the other hand, the rotation axis is inclined 
from the magnetic field lines at an angle of Sq — 60° in 
models 3 and 4. Model 3 has a weaker magnetic field 
than model 2, whereas model 4 has the strongest mag- 
netic field among all the models, as shown in Table 12.11 
In addition, the rotation rate of model 4 is slightly slower 
than those of models 2 and 3. 

4.1. Magnetohydrodynamic Simulations 

In this subsection, we simply show the results of MHD 
simulations. As described in §]2T] we calculated the 
evolution of the magnetized collapsing cloud from the 
prestellar cloud stage. Figure [3] shows the cloud evolu- 
tion for model 1, in which the density distribution (colors 
and contours) and velocity vectors on the y = cutting 
plane of the I = 7 grid are plotted. Figure Eta) shows the 
initial state. Figure |3jb) shows the collapsing cloud be- 
fore protostar formation (i.e., the prestellar stage), and 
Figures [3{c) and (d) show the cloud after protostar for- 
mation (i.e., the protostellar stage). In this model, a 
protostar forms at t = 7.1 x 10 4 yr after the cloud begins 
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Model 1 : No rotation 




magnetic field lines 



Model 3 : 

Non-aligned rotation, 
weak B 




Model 2 ; Aligned rotation 

rotational axis 



Model 
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Non-aligned rotation, 
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magnetic field lines 



Fig. 3. — Density distribution (colors and contours) and velocity 
vectors (arrows) on the y = cutting plane at four different epochs 
for model 1. Elapsed time t and protostellar mass M ps are listed 
in each panel. 



Fig. 2. — Schematic view of the initial conditions for each model. 
Blue and red arrows represent the directions of initial magnetic field 
lines and initial rotation axis, respectively. The width of the blue 
arrows reflects the magnetic field strength. The magnetic field lines 
are aligned with the rotation axis for model 2. whereas they are 
inclined from the rotation axis by <5n = 60° for models 3 and 4. 

to collapse. At the epoch corresponding to Figure [HJd) , 
the protostar has a mass of 0.55 M and is still embed- 
ded in the collapsing cloud. Because the initial cloud in 
model 1 has a magnetic field but no rotation, the gas col- 
lapses essentially along the magnetic field lines. Thus, a 
disk-like structure (pseudo-disk) forms, extending in the 
direction perpendicular to the magnetic field lines, as 
seen in Figure E)Jd). 

Figure 2] shows snapshots of the protostellar stage for 
models 2, 3, and 4. After protostar formation, model 

2 shows a weak outflow that is traced as a shock front 
near z ~ 800 AU, far from the center of the cloud in Fig- 
ure HJa). Note that we cannot resolve the outflow driven 
by the circumstcllar disk with sufficient spatial resolu- 
tion because we adopted sink cells in the high-density gas 
region in order to perform a long-term calculation (see 
£|5.4[) . The large- and small-scale structures for model 

3 are plotted in Figures |U[b) and (d). For this model, 
the disk is not perpendicular to the z-axis, that is, the 
disk normal is inclined from the direction of the global 
magnetic field lines. This is because the rotation axis 
of the initial cloud is inclined from the initial magnetic 
field lines. In addition, for model 3, the initial rotational 
energy (/3o = 0.02) is equivalent to the initial magnetic 
energy (70 = 0.02). Thus, the cloud's rotation can af- 
fect the dynamical evolution of the cloud and the disk 
is inclined from the z-axis, as seen in Figure 0Jd). On 
the other hand, although the rotation axis is inclined 
from the magnetic field direction in the initial cloud, the 



disk tends to form in the direction perpendicular to the 
magnetic field lines for model 4 (Figure [4jc] ) . For this 
model, because the initial magnetic energy (70 = 0.57) 
dominates the initial rotational energy (/3n = 0.01), the 
cloud rotation is expected to have very lit tle effect on the 
dyna mical cloud evolution (for details, see Machida ct al. 
|2006j). In addition, because of the strong magnetic field, 
magnetic braking is more effective in model 4 than in 
models 2 and 3. 

For all the models, over half of the initial cloud mass 
falls onto either the protostar or the circumstcllar disk 
by the end of the calculation. This means that we calcu- 
lated the cloud (or protostar) evolution until the Class 
or I protostar stage. We will show the evolution of the 
magnetic field configuration in the next subsection. 

4.2. Expected Distributions of Polarization 

After MHD simulations were performed, we calculated 
the polarization of the thermal dust e mission accordin g 
to the procedure described in £l3l and rTbmisakal (|2011h . 
In this subsection, we show the expected distributions 
of the polarization. Some of the arguments in the next 
subsection are simila r to p revious studies (|Frau et al.l 
l201lHGoncalves et al.|[2008l) . However, to compare mod- 
els with and without rotation, or compare aligned rota- 
tion model with misaligned rotation model, we describe 
polarization patterns for all models in detail. 

4.2.1. Model 1: Non-rotating Cloud 

Figure [5] plots the polarization distribution for model 
1 at different epochs viewed from the angle of (0,4>) = 
(90°, 0°). In this model, the initial cloud is assumed 
to have a uniform magnetic field but no rotation. In 
the figure, the colors, black contours, and straight lines 
represent the polarization degrees, column densities, and 
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FlG. 4. — Density distribution (colors and contours) and velocity 
vectors (arrows) on the y = cutting plane for models 2 (a), 3 (b 
and c), and 4 (c). Elapsed time t and protostellar mass Mp S are 
listed in each panel. Panel (d) is a close-up view of panel (b). 



0=90* 
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-MOO -2000 -1000 1000 5000 3000 -3000 -2000 -1000 1000 2000 3000 

am 5(au) 



Fig. 5. — Polarization maps for model 1 viewed from the angle 
(9, <f>) = (90°, 0°) at four different epochs, which correspond to the 
respective panels of Figure|3] The polarization degree (colors), po- 
larization vector (lines), and column density (contours) are plotted 
in each panel. The number density at the center of the cloud n c 
or protostellar mass M pB is shown in each panel. 



polarization vectors (the direction of the B-vector of the 
electromagnetic waves), respectively. Hereafter, we call 
this type of figure a polarization map. Viewed from 
(9, 4>) = (90°, 0°), the £ and r\ directions in the obser- 
vation grid correspond to the y and z directions in the 
simulation grid, respectively (see Figure [lj. Thus, each 
panel of Figure [5] corresponds to the same panel in Fig- 
ure [3] (Note that the distributions on the y = and 
x = planes are the same because the system is axisym- 
mctric. In addition, because this model maintains the 
axisymmctry, the physical quantities do not depend on 
the azimuthal viewing angle 0. Thus, we show only the 
polarization map with = 0° for this model.) 

Figure Eta) is a polarization map of the initial state for 
model 1, in which all the polarization vectors are aligned 
with the 77-axis. Figure EJb) shows the polarization just 
before protostar formation. At this epoch, the polariza- 
tion vectors are distorted near the center of the cloud, 
whereas the column density contours maintain a spher- 
ical symmetry Figure EJc) shows the polarization after 
protostar formation. The polarization degree around the 
protostar (or around the center of the cloud) gradually 
decreases because the magnetic field parallel to the line 
of sight begins to contribute as the cloud collapses. Af- 
ter protostar formation, the magnetic field lines point to 
the protostar around the center of the cloud, generating a 
magnetic field parallel to the line of sight. As the parallel 
component increases, the polarization degree decreases, 
because only the magnetic field lines perpendicular to the 
line of sight contribute to the polarization. (We assume 
that the magnetic field changes its direction while main- 
taining its strength.) Around the center of the cloud, 
the polarization degree decreases to ~ 13% when the 
protostellar mass is M ps = 0.05 M & (Figure [5[c] ) , and 
to < 10% when M ps = 0.5 M@ (Figure EM). This indi- 



cates that the magnetic field lines are gradually inclined 
toward the line-of-sight direction. 

From the polarization vectors in the £77 plane, we can 
directly confirm that the magnetic field lines are inclined 
from their initial direction. In other words, the magnetic 
field lines point to the protostar around the center of 
the cloud. Thus, we can imagine an hourglass structure 
of the magnetic field lines on the polarization maps. In 
addition, the column density contours show a disk-like 
configuration, which indicates the existence of a pseudo- 
disk, as described in £14.11 

To investigate the effects of viewing angle 9 on the 
polarization, we selected a snapshot when the protostar 
has a mass of M ps = 0.5 M@. Figure El shows polariza- 
tion maps with different viewing angles of 9 = 0° (a) , 30° 
(b), 45° (c), 60° (d), 80° (e), and 90° (f); the azimuthal 
angle is fixed at = 0°. The viewing angle 9 = 0° cor- 
responds to the pole-on view and 9 = 90° corresponds 
to the edge-on view. Thus, Figure Etf) is identical to 
Figure EJd) . Figure El indicates that the properties of 
the polarization map depend strongly on the viewing an- 
gle 9. The figure shows that the polarization degree of 
the 9 = 0° map (panel [a]) is considerably lower than 
those for 9 7^ 0°. We assumed that the initial mag- 
netic field lines arc parallel to the z-axis. Thus, as the 
cloud evolves, the z component of the magnetic field (B z ) 
dominates the other components (B r and B^) in the en- 
tire cloud. Therefore, the lower polarization is a natural 
consequence of the configuration when we measure the 
polarization from the pole (or we integrate the polariza- 
tion along the z-axis) because the magnetic field lines are 
almost parallel to the line of sight. In Figure Et a )> the 
central region has a larger polarization degree, which in- 
dicates that magnetic field lines perpendicular to the line 
of sight exist around the center of the collapsing cloud. 
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In MHD simulations, the azimuthal component of the 
magnetic field (B^) is not generated because model 1 
has no rotation, whereas the radial component (B r ) is 
generated as the cloud collapses. Because the magnetic 
field lines are strongly distorted near the protostar, the 
central region has a larger polarization degree owing to 
the non- negligible component of B r in Figure EJa) . 

On the polarization maps for 9 = 30°, 45°, and 60° 
in Figures EJb), (c), and (d), respectively, an interesting 
feature appears around the protostar. In these maps, a 
region of very low polarization with a polarization degree 
of < 2% (white region) appears in the transverse direc- 
tion around the protostar (or along the £-axis). The low- 
polarization area shrinks as the viewing angle 9 increases. 
Finally, the region disappears on the polarization maps 
of 6 = 80° (Figure [le]) and 90° (Figure [Bp]). Thus, 
this low polarization is attributed to the projection ef- 
fect of the magnetic field lines. To compare the projected 
magnetic field lines on an arbitrary plane with those in 
three dimensions, a three-dimensional view of the mag- 
netic field lines is plotted in Figure [7] Figure [3a) shows 
the magnetic field lines before protostar formation (cor- 
responding to Figure [SJb] ) , whereas panels (c) and (d) 
show those after protostar formation (corresponding to 
Figure [SJd] ) . In addition, panel (d) has a viewing angle 
of 6 = 60°, whereas the left panels, (a) and (c), have 
9 = 90° (i.e., an edge-on view). 

The polarization degree is lower around the protostar 
because of two reasons. First, the magnetic field parallel 
to the line of sight does not contribute to the integration 
of the polarization. The magnetic field lines are strongly 
bent toward the protostar around the center of the cloud. 
As seen in Figure [3 these field lines have a component 
parallel to the line of sight. Thus, the polarization de- 
gree decreases near the center of the cloud. The other 
reason is the cancelation of positive Stokes (q, u) and 
negative Stokes (<?, u) along one line of sight (see Fig- 
ure E[c]). That is, two overlapping magnetic field lines 
cancel the polarization if they are perpendicular to each 
other because the direction of dust alignment is perpen- 
dicular. The model has an axisymmetric magnetic field 
with respect to the z-axis and a line-symmetric distri- 
bution with respect to the z = plane. In this case, 
as shown in Figure EJd), the foreground magnetic field 
(red) and the background field (blue) overlap each other 
with an angle of <; 90°. This weakens the polarization in 
the polarization map. This occurs in the region z ~ or 
j] ~ in a symmetric way with respect to the z-axis or 
7/-axis in this case. 

In the collapsing cloud, the magnetic field lines have 
an hourglass-shaped configuration. As shown in Fig- 
ures [7|a) and (c), a clear hourglass appears when we 
see the magnetic field lines from the edge (the edge-on 
view, 8 = 90°). On the other hand, a complicated struc- 
ture is expected around the center of the cloud when 
viewed from 9 < 90° (Figure [7[d] ) . In reality, we can- 
not observe the three-dimensional structure of the mag- 
netic field lines; we observe only the projected magnetic 
field onto the polarization map. Instead, Figure [6j indi- 
cates that we can roughly estimate the angle between the 
global magnetic field lines and the line of sight from the 
distribution of the low-polarization region. 

We also investigated the scale dependency of the con- 



figuration of the magnetic field lines. Figure [5] shows po- 
larization maps with 9 = 90° at different spatial scales. 
Panel (a) is 1 2000 AUxl 2000 AU in size, whereas panel 
(b) is 3000 AUx 3000 AU in size. Figure M clearly indi- 
cates that the polarization vectors closer to the protostar 
are more distorted. Thus, the hourglass structure is em- 
phasized at smaller scales. 

4.2.2. Model 2: Rotating Cloud 

To investigate the effect of cloud rotation on the po- 
larization, we calculated the polarization for model 2, 
which has the same magnetic field strength as model 1 
but a finite angular momentum. In model 2, the rota- 
tion axis is parallel to the magnetic field lines, as shown 
in Figure [2] Figure [9] shows the polarization maps for 
model 2 at different epochs for (0, 0) = (90°, 0°). Panels 
(a) and (b) show the polarization before and just after 
protostar formation, respectively. At these epochs, the 
polarization distribution for model 2 is almost the same 
as that for model 1 (compare Figure[9]with Figure[5|). All 
the panels of this figure indicate an hourglass configura- 
tion of the magnetic field lines. In addition, a region of 
slightly low polarization (orange region) appears above 
and below the center of the cloud, as seen in Figure [5j 
Thus, the cloud rotation has very little effect on the po- 
larization distribution for the edge-on view before and 
just after protostar formation. 

On the other hand, after protostar formation, the po- 
larization degree of model 2 is very different from that of 
model I with no rotation. As shown in Figure [9jc), in the 
rotating cloud, the polarization degree just above and be- 
low the protostar is extremely low (< 2%, white region). 
Even in model I (non- rotating cloud), a low-polarization 
region appears around the protostar (Figures [5jc] and 
[d]). However, the polarization degree around the pro- 
tostar is ~ 10% in model 1 but < 2% in model 2. In 
the non-rotating model (model 1), the low polarization 
degree is caused by the generation of the radial com- 
ponent of the magnetic field, which produces magnetic 
field lines parallel to the line of sight when viewed edge- 
on. On the other hand, in the rotating cloud, in addition 
to the radial component, a toroidal component (B^,) is 
generated by the rotation. Because the system is axisym- 
metric with respect to the z-axis, the angle between the 
foreground and background magnetic field lines is large 
(~ 90°) when the magnetic field consists of poloidal and 
toroidal components of comparable strength (see Fig- 
ure E[d]). Consequently, a region of extremely low po- 
larization appears above and below the protostar, where 
the toroidal component is amplified to become compa- 
rable with the poloidal one. This region appears only 
after protostar formation because the toroidal magnetic 
field is a result of rotational motion around the proto- 
star. Note that in the isothermal collapse phase (i.e., 
before protostar formation), such a strong toroidal field 
is not expected because the collapse timescale is shorter 
than the rotational timescale. Therefore, the extremely 
low polarization degree around the center of the cloud is 
an indicator of the existence of the protostar and rotating 
disk. 

Figure [TU] shows polarization maps with different view- 
ing angles (9 = 0° [a], 30° [b], 45° [c], 60° [d], 80° [e], 
and 90° [f]; = 0°) for model 2 when the protostar has 
a mass of M ps = 0.5 M Q . Similar to model 1, because 
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Fig. 6. — Polarization maps for model 1 viewed from 8 = 0° (polc- 
f); </> = 0° when M ps = 0.5 Mq. 



(a) t-7.0 X 10 J [yr] 6-90° (b) depolarization mechanism 




Fig. 7. — Magnetic field lines in three dimensions before (a) and 
after (c) protostar formation for model 1. In panel (d), 8 = 60°, 
and the two panels on the left show the edge-on view (8 = 0°). 
Red lines represent magnetic field lines traced from the lower and 
upper boundary surfaces. White contour lines and colors repre- 
sent the density distribution on the y = cutting plane. Panel 
(b) gives a schematic explanation of the depolarization mechanism 
when we observe the cloud from 8 = 60°. In the region marked 
in panel (d), the foreground and background magnetic field lines 
intersect at a large angle ~ 90°. In this case, dust alignment be- 
tween the foreground and background is canceled, which decreases 
the polarization degree. 

model 2 maintains axisymmctry with respect to the z- 
axis, the polarization distribution does not depend on 
the azimuthal viewing angle <f>. A comparison of Fig- 
ures [TU] (rotating model) and [S] (non-rotating model) re- 
veals a difference in the polarization distribution. In the 
non-rotating model, a low-polarization region (< 2%) ap- 
pears in a transverse direction (^-direction) around the 



view; a), 30° (b), 45° (c), 60° (d), 80° (e), and 90° (edge-on view; 
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Fig. 8. — Polarization maps for model 1 at different spatial scales 
with (0, <f>) = (90°, 0°) when M ps = 0.5 Mq . Right panel shows 
the central region in the left panel enlarged four times. 



protostar at viewing angles of 9 = 30°, 45°, and 60° 
(Figure H3) . This feature was explained using Figure [7] in 
i )4.2.1l On the other hand, in the rotating model, the low 
polarization degree appears in two regions, upper-right 
and lower-left regions of the center, as shown in the po- 
larization maps with 9 = 30°, 45°, and 60° in Figure [TU1 
Thus, in model 2, the line symmetry of the polarization 
distribution with respect to the 77-axis is broken and a 
point symmetry with respect to the protostar appears. 
This point-symmetric distribution of t he polarization i s 
a striking featu re of rotating clouds ()Tomisakal 120111) . 
iTomisakal (|2011h pointed out that (i) in the protostellar 
phase, the toroidal magnetic field arises owing to rota- 
tion of the central disk. This toroidal component has 
antisymmetry with respect to the z = plane (for exam- 
ple, B,p > for z > 0, and B$ < for z < 0). (ii) The 
magnetic field consisting of the toroidal and hourglass- 
type poloidal components becomes point-symmetric (see 
his Figure 10). When a magnetic field of this type is seen 
from 30° <, 9 <^ 60°, the region with a low polarization 
degree is distributed in a point-symmetric way. In other 
words, this is caused by a combination of viewing angle 
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Ca) t=7.2 x 10'[yr]. n-3.8 x 1 B [cm '] 

,77 



(b) t=7.5x ipy]. M n ^0.01[Msun] 




Fig. 9. — Polarization maps for model 2 at three different epochs 
(t = 7.2 X 10 4 yr, 7.5 X 10 4 yr, and 1 X 10 5 yr) at the same viewing 
angle of (6>, <f>) = (90°, 0°). Panel (a) represents the prestellar 
stage, and panels (b) and (c) show the protostellar stage. The 
protostar masses at each snapshot are M pB = 0.01 Mq (panel [b]) 



and M„ 



0.5 Mq (panel [c]). 



and rotation effects. 

As described above, a strong toroidal field also causes 
a low polarization degree. This low polarization appears 
in the region just above and below the protostar and 
exhibits line symmetry with respect to the £-axis, as 
shown in Figure fTOl f). Thus, a low-polarization region 
just above and below the central protostar appears in 
the 9 = 45°, 60°, 80°, and 90° maps in Figures [TO^c), 
(d), (e), and (f), respectively, and is caused by the coex- 
isting toroidal and poloidal magnetic field components. 
In summary, two effects decrease the polarization de- 
gree: line-symmetric and point-symmetric effects. Both 
arise from cancelation between the foreground and back- 
ground dust alignments. Even without the toroidal field 
(model I), cancelation in the poloidal component pro- 
duces a low-polarization region extending near the £-axis. 
On the other hand, a magnetic field consisting of poloidal 
and toroidal components causes vertical extension of the 
low-polarization region near the 77-axis. In this case, the 
former cancelation occurs in the upper-right and lower- 
left regions of the center in a point-symmetric way when 
viewed at angles in the range of 30° <^ 9 <^ 60°. 

To clarify the effects described above, we plot a three- 
dimensional view of the magnetic field lines in Fig- 
ures ITTT a) and (b) at viewing angles of = 90° and 
9 = 60°, respectively. The figure clearly shows that the 
toroidal field is confined in a small region around the pro- 
tostar, and the hourglass configuration of the magnetic 
field lines extends to the large scale. 

Although the polarization distribution in the rotating 
model differs from that in the non-rotating model (com- 
pare Figure [6] with Figure [10]) , the polarization vectors 
and column density contours of the models are almost 
the same at the large scale. In addition, the hourglass 
configuration at the cloud scale can be confirmed from 
the polarization maps in both models. At the small 



scale, however, the polarization vectors of the rotating 
model deviate from the hourglass configuration near the 
protostar because the toroidal component dominates the 
poloidal one. In terms of observations, it is difficult to 
observe the deviation because the polarization degree is 
too low to allow identification of the direction of the po- 
larization in such a region. 

4.2.3. Model 3: Misaligned Cloud with Weak Magnetic 

Field 



In §4.2.2[ we investigated the polarization distribution 
in a cloud with an idealized initial setting in which the 
rotation axis is parallel to the magnetic field lines. In re- 
ality, however, the rotation axis is not necessarily parallel 
to the magnetic field lines in molecular cloud cores. In 
this subsection, we investigate the polarization in model 
3 in which the initial cloud has a rotation axis that is 
not parallel to the magnetic field lines (sec Figure [2]) . In 
model 3, the initial angle between the rotation axis and 
the magnetic field lines is So = 60°, and the rotational 
energy (/?o = 0.02) is equivalent to the magnetic energy 
(70 = 0.02) in the initial cloud. Figure [T2] shows po- 
larization maps at two viewing angles, (0, 4>) = (90°, 0°) 
([a], [b], and [c]) and (0,0) = (90°, 90°) ([d], [e], and [f]), 
for model 3 at three different epochs, in which two snap- 
shots (t = 6.3 x 10 4 yr and t = 7.3 x I0 4 yr) describe the 
prestellar phase ([a], [b], [d], and [e]) and one snapshot 
(t = 1.1 x 10 5 yr) represents the protostellar phase ([c] 
and [f]). At an age of t = 1.1 x 10 5 yr, the protostar has 
grown to M ps = 0.5 Mq. As in equation ©, the rota- 
tion axis is inclined by So = 60° from the z-axis on the 
xz plane in the simulation grid. Because the yz plane 
in the simulation grid coincides with the £77 plane in the 
observation grid when viewed from (9, (j>) = (90°, 0°), the 
initial rotation axis is inclined by So = 60° from the 77- 
axis to the observer in the upper panels of Figure [T2l (sec 
Figure[l|). In the lower panels, where (0, 0) = (90°, 90°), 
the initial angular momentum vector is directed toward 
(C,'7) = (-v / 3/2,+l/2). 

Figure [T^] shows that the polarization degree around 
the center of the cloud decreases with time because the 
low-polarization region develops after the protostellar 
grows sufficiently ([c] and [f], in which M ps = 0.5 Mq). 
The reduction rate in the polarization degree for model 
3 (Figure [T2|) is similar to that for model 2 (Figure |9]) . 
However, a line symmetry with respect to the 77-axis is 
broken in model 3 (panel [b]) but maintained in model 
2 (Figure ^ . In addition, the polarization vectors in 
Figures IT2T b) and (c) differ considerably from those in 
Figures [5] and [9] For models 1 and 2, the magnetic field 
lines have an hourglass configuration for any viewing an- 
gle </> (and 0), indicating that they have axisymmctry 
around the z-axis. On the other hand, for model 3, the 
polarization vectors appear different at viewing angles 
of <j> = 0° (b) and <j> = 90° (e). Panels (a) and (d), 
which represent the structure before core formation, in- 
dicate a slight asymmetry with respect to the 77-axis at 
both viewing angles. This means that in the prestellar 
phase, the departure from axisymmetry around the z- 
axis is weak. Panels (c) and (f), which show polarization 
maps of the protostar phase with M ps = 0.5 Mq, seem to 
indicate that the axisymmetry in the polarization vectors 
is recovered again. Therefore, the departure from the ax- 
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Fig. 10. — Polarization maps for model 2 from the viewing angles of 9 = 0° (pole-on view; a), 30° (b), 45° (c), 60° (d), 80° (e), and 90° 
(edge-on view; f); = 0° when M ps = 0.5 Mq. 
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Fig. 11. — Magnetic field lines in three dimensions when M pB = 
0.5 Mq for model 2. Red lines represent magnetic field lines. White 
contours and colors represent density distribution on the y = 
cutting plane. The viewing angles are 9 = 90° (a) and 9 = 60° 
(b). 

isymmetric structure is prominent in the final phase of 
prestellar contraction (panels [b] and [c]) and the early 
protostellar phase (Figure I13[) . This non-axisymmetric 
magnetic field has the shape of the letter S in panel (b) . 
That is, in panel (b), the polarization vectors change 
their directions from the vertical direction in the upper- 
most region to the direction from the upper left to the 
lower right around the equator (j) = 0), and finally to 
the vertical direction in the lowermost region. The non- 
axisymmetric magnetic field appears in another way in 
panel (c): the major axis of the total intensity distribu- 
tion (column density) extends from the upper right to 
the lower left. Although the polarization vector is essen- 
tially perpendicular to the disk major axis in models 1 
and 2, the polarization vector in panel (e) indicates an 
hourglass shape, but its axis is not perpendicular to the 
major axis of the disk. 

As described in ^4.2. 2[ the polarization vectors do not 
depend on the azimuthal viewing angle in a cloud with 
an initial rotation axis parallel to magnetic field lines. 
On the other hand, they do depend on <p when the ini- 



tial magnetic field lines are misaligned with the rotation 
axis (Figure [12]) . Figure [13] shows polarization maps for 
9 = 90° for model 3 at the epoch when the protostar has 
a mass of M ps = 0.1 Mq and a three-dimensional view 
of the magnetic field and isodensity surface (panel [b]). 
The viewing angle <j> is equal to 0° in panel (a) and 90° 
in panel (c). The protostellar mass of M vs = O.1M 
is smaller than M ps = O.5M in Figures 1127 c) and (f). 
This clearly shows that in the early evolutionary stage 
of protostar accretion, the viewing angle has a signifi- 
cant effect. Note that the early protostellar phase (Fig- 
ures [T3Ia] and [c]) is very similar to the final phase of the 
prestellar stage (Figures [T2]b] and [e]). The polarization 
vectors have the S'-shaped configuration in Figure [T3Ta) 
at 4> = 0°, whereas they have an hourglass configura- 
tion in Figure [T3T c) at <fi = 90°. Although the disk is 
inclined with respect to the z-axis (the isodensity sur- 
face in Figure I13|b] and the surface density contours in 
Figure I13[c] ) , the axis of the hourglass-shaped magnetic 
field seems parallel to the z-axis. In addition, each panel 
shows a different distribution of the polarization degree. 
A low-polarization region appears in the lower-right and 
upper-left regions toward the protostar for = 0° (panel 
[a]), whereas it appears almost along the £-axis near the 
protostar for <j> = 90° (panel [c]). In reality, for model 
3, the magnetic field lines do not show a clear hourglass 
structure in three dimensions. Nevertheless, depending 
on the viewing angle, they appear as an hourglass shape, 
as shown in Figure 113( c) . 

The distribution of the polarization degree and polar- 
ization vectors for <fi = 90° seems to be almost the same 
as that in the aligned rotation case (model 2). However, 
the column density contour in Figure H~3T c) indicates that 
a pseudo-disk is inclined by ~ 45° from the £-axis. This is 
an outcome of the inclined density distribution shown in 
Figure [TST b). That is, when viewed from the y-axis, the 
disk (isodensity surface) also extends from the lower left 
to the upper right. On the other hand, the column den- 
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Fig. 12. — Polarization maps for model 3 at different epochs viewed from (9, <f>) = (90°, 0°) (a, b, and c) and (0, <j>) = (90° ,90°) (d, e, 
and f). Panels (a) and (d) (f = 6.3 X 10 4 yr) and panels (b) and (e) (t = 7.3 X 10 4 yr) correspond to the prestellar contraction phase. In 
contrast, panels (c) and (f) cover the protostcllar phase with a protostar mass of 0.5 Mq (t = 1.1 X 10 5 yr). 
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Fig. 13. — Polarization maps when M ps = 0.1 (a and 
c) for model 3. Viewing angles arc (9,tf>) = (90°, 0°) (a) and 
(9,<t>) = (90°, 90°) (c). Panel (b) represents the three-dimensional 
distribution of the magnetic field lines and isodensity surface of 
n = 4.7 X 10 5 cm~ 3 . 

sity has a spherical symmetry in Figure [T5ta) for = 0°. 
Figure fT5Tb) shows that the symmetry rotation axis of 
the disk is in the xz plane, and the disk normal is directed 
to approximately ~ 45° from the x-axis. 

The deviation from the hourglass configuration occurs 
because the angular momentum is not parallel to the 
global magnetic held. However, as the cloud collapses, 
the angular momentum perpendicular to the global mag- 
netic fiel d lines is effectivel y transferred by magnetic 
braking (|Price fc Bate! 120071 ). Thus, at a later evolu- 
tionary stage, the rotation axis tends to align with the 



global magnetic field lines. Therefore, the rotation axis 
becomes parallel to the magnetic field, and the hourglass 
configuration is recovered. However, because the angular 
momentum of the rotating disk is not completely trans- 
ferred, the disk normal is not completely aligned with the 
global magnetic field, especially at the small scale (~ 500 
AU) , as seen in the central regions of Figures [T2T c) and 
(0- 

Cloud rotation has a significant effect on the star for- 
mation process because it produces a circumstellar disk 
where planets form. It is thought that magnetic field 
lines are generally not well aligned with the rotation axis 
in molecular cloud cores. Thus, we expect that the S- 
shaped configuration of the magnetic field lines seen in 
Figuresfr^land [T3]will be frequently observed in the early 
stage of star formation in future high-angular-rcsolution 
observations. 

4.2.4. Model J t : Misaligned Cloud with Strong Magnetic 

Field 

As described in £|4.2.3[ when the rotation axis is in- 
clined from the magnetic field lines, the polarization 
vectors do not indicate hourglass-shaped but S-shapcd 
structure. The S'-shaped configuration is caused by cloud 
rotation, which can cause the magnetic field lines to devi- 
ate from the hourglass configuration. Both the magnetic 
field (Lorentz force) and rotation (centrifugal force) can 
form a disk in the collapsing cloud. When a gas cloud 
having a strong magnetic field collapses along the mag- 
netic field lines, a pseudo-disk forms in the direction per- 
pendicular to the magnetic field lines. On the other hand, 
when a gas cloud having a larger angular momentum col- 
lapses along the rotation axis, a rotating disk is extended 
in the direction perpendicular to the rotation axis. Thus, 
it is expected that the former creates a clear hourglass 
polarization configuration and the latter shows the S'- 
shaped configuration, depending on the viewing angle, 
when the magnetic field lines are not parallel to the ro- 
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Fig. 14. — Polarization maps for model 4 when M pB = 0.01 Mq 
at viewing angles of (0, <j>) = (90°, 0°) (a) and (9,tf>) = (90°, 90°) 
(b). 

tation axis. Wc confirmed the S'-shaped configuration 
with a relatively weak magnetic field in i )4.2.3l (i.e., the 
former case). In this subsection, to investigate the lat- 
ter case, we show polarization maps of model 4. Model 4 
initially has a strong magnetic field (70 = 0.57) but weak 
rotation (/3 = 0.01). 

Figure [14] shows polarization maps of model 4 viewed 
from the angles of (0,0) = (90°, 0°) (a) and (0, </>) = 
(90°, 90°) (b) just after the protostar formation epoch 
(M ps = 0.01 Mq). The figure indicates that polarization 
vectors have the hourglass configuration irrespective of 
the viewing angle <j). On the other hand, in model 3, the 
polarization vectors exhibit the S-shaped configuration 
even when M ps = 0.1 Mq. We confirmed that for model 
4, the hourglass configuration is maintained until a later 
evolutionary stage. Thus, we conclude that in a cloud 
with a relatively strong magnetic field, rotation rarely 
causes the magnetic field lines to deviate from the hour- 
glass configuration. This is because gas falls onto the 
central region mainly along the magnetic field lines, not 
along the rotation axis. 

At a small scale, however, the line symmetry of the 
polarization vectors with respect to the £-axis seems to 
be broken. In addition, the column density contours in 
Figure [14] imply an inclined disk. Thus, the effect of 
the misalignment between the magnetic field and the 
rotation axis cannot be negligible even for this model. 
The deviation of the magnetic field lines from the hour- 
glass configuration very close to the protostar may be 
detectable with future instruments such as ALMA. 

5. DISCUSSION 

5.1. Polarization Angle and Configuration of Magnetic 
Field Lines 

To date, magnetic field lines are generally considered 
to have an hourglass structure in a gravitationally con- 
tracting cloud. However, as described in £|4.2[ we showed 
that this hourglass configuration is not always realized 
in a collapsing cloud. Depending on the viewing angle, 
the magnetic field lines have an S-shaped configuration 
when the cloud rotation influences the dynamical cloud 
evolution. In any case (hourglass or S'-shaped configu- 
rations), the magnetic field lines gradually deviate from 
their initial configuration as the cloud collapses. In this 
subsection, wc present a method of measuring the mag- 
netic field structure in star-forming clouds by qualita- 
tively estimating the magnetic field line structure. 

In the low-density gas region (n < 10 u cm- 3 ), the 
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Fig. 15. — Time sequence of deviation angle of polarization vec- 
tors from the initial configuration for models 1 (a and b) and 3 
(c and d). Deviation angle was calculated along the £-axis at 
r\ = 253 AU for (a) and (c) and at 77 = 860 AU for (b) and (d) 
on the polarization map with (8, <j>) = (90°, 0°). Viewing angle 
dependence is shown in panels (e) and (f), in which snapshots of 
model 3 at t = 7.7 X 10 4 yr are shown for <f> = 0°, 30°, 60°, and 
90°. 



magnetic field is well coupled with neutrals via momen- 
tum exchange with ions. Thus, as the cloud collapses, 
the magnetic field lines converge toward the center of 
the cloud and form the hourglass structure, when they 
are initially parallel to the rotation axis. If we simply as- 
sume an initially uniform magnetic field, the transverse 
motion of the gas drags the magnetic field lines. In a 
cloud, the gas initially distributed at a distance from the 
cloud center must move a long distance and later fall onto 
the cloud center (or onto the protostar or circumstellar 
disk). Thus, the deviation of the magnetic field from the 
initial configuration is believed to propagate from the 
center to the outer envelope with time. To qualitatively 
investigate this, we calculated the deviation angle of the 
polarization vectors from the initial configuration on the 
polarization map with (0» = (90°, 0°) 6 . Figure [TBI 
shows a time sequence of the deviation angle in which 
the horizontal axis coincides with the £-axis. The angle 
is set to 0° for the initial magnetic field lines. The devi- 

6 We owe the idea of figure [TBI to S.Basu, who presented similar 
diagrams in Winter School and Worksh op on Star Formation in 
Tokyo, 2011 (see also lBasu eTaLl I l2009al lbr0 



Magnetic Field Structure in Star-Forming Cores 



13 



ation angle is measured along a straight line parallel to 
the £-axis at r\ = 253 AU in panels (a), (c), and (e) and 
860 AU in panels (b), (d), and (f). 

Figures HSJa) and (b) show the deviation angle for 
model 1. Model 1 has no rotation and thus exhibits a 
clear hourglass, as shown in Figure [5] The snapshots are 
the same in Figure \S\ and Figures ITSTa) and (b). Fig- 
ure I15f a) shows that the deviation angle has two sharp 
extremes at r ~ ±(100 — 200) AU. The maximum devia- 
tion angle is ~ 10° before protostar formation and <J 40° 
after it. That is, the deviation angle gradually increases 
overall. This indicates that the magnetic held lines are 
gradually inclined from their initial configuration. In ad- 
dition, the deviation angle has a point symmetry with 
respect to the origin. The point symmetry in this plot is 
an indicator of the hourglass configuration of the mag- 
netic field lines. Figure ITST a) shows that the deviation 
angle at a larger distance (> 1000 — 2000 AU) remains 
< 5°. This is because the deviation angle in this panel is 
measured along the line near the equator. As shown in 
Figure [5J the polarization vectors arc inclined very little 
near the r/ = symmetry axis; this is because the mag- 
netic field lines on the equatorial plane are pulled mainly 
in the direction perpendicular to their initial direction 
(or in the direction toward the protostar). 

A comparison of panels (a) and (b) indicates that the 
absolute values of the extremes decrease as the height 
from the mid-plane increases. The reason is that the 
magnetic field is dragged by the gas, which is contracting 
toward the center. However, at a distance from the center 
in the £ direction, the deviation angle is larger for a larger 
height (panel [b] ) . This seems to be explained by the fact 
that because of the mirror symmetry with respect to the 
z = plane, the magnetic field near the z — plane must 
be aligned along the z-axis. Away from the z = plane, 
the magnetic field is bent more strongly. 

Figures I15f c) and (d) , in which the snapshots are the 
same as in Figures [12] and I13[ show that the point sym- 
metry is broken in model 3. In model 3, the initial 
rotation axis is inclined from the global magnetic field 
lines, and the .S'-shaped configuration of the polariza- 
tion vectors appears (t = 7.3 x 10 4 yr in Figure I12|b] 
and t = 7.7 x 10 4 yr in Figur <JT3|) . In the panel, a single 
negative extreme appears at £ ~ 200 — 300 AU in the 
early evolutionary stage, which represents the S'-shaped 
configuration of the polarization vectors (Figure I12|b] ) . 
On the other hand, another positive extreme develops 
gradually with time, and the deviation angle begins to 
have two extremes (one positive and one negative) . This 
means that the magnetic field configuration changes from 
the S-shaped to the hourglass with time. This deviation 
angle plot can clearly distinguish between the hourglass 
and S-shaped configurations. 

Figures IT5T e) and (f) show the viewing angle depen- 
dence of 4>. Viewed from = 90°, the deviation angle of 
model 3 has two extremes, one positive and another neg- 
ative. However, the same model shows a single extreme 
if we observe from <fi < 60°. This clearly shows that 
in these periods the S'-shaped configuration (<f) <, 60°) 
and the hourglass configuration (<fi ~ 90°) arc realized 
simultaneously and that the appearance depends on the 
viewing angle <p. 



5.2. Toroidal Components on Polarization Map 

In a non-rotating collapsing cloud, only the poloidal 
component of the magnetic field is generated. In con- 
trast, the toroidal field appears in addition to the 
poloidal field in rotating clouds. The protostellar outflow 
and angular momentum transfer due to magnetic brak- 
ing ar e closely related to the cloud rotation and toroidal 
field (|Tomisakal 120021 : iMachida et all 1201 lh . Thus, the 
toroidal field is crucial to identify the driving mecha- 
nism of protostellar outflows or investigate the transfer 
mechanism of angular momentum. In this subsection, 
we present the observation of toroidal components on 
the polarization maps. 

When we observe the star-forming cloud from the pole 
(pole-on view) , we can directly observe the toroidal field. 
We can confirm the presence of the toroidal component 
of the polarization vectors in Figure [TUta) (pole-on view 
of polarization map). However, the polarization degree is 
very low at the large (or cloud) scale because the poloidal 
component dominates the toroidal component. Note that 
because the frecfall timcscale is shorter than the rotation 
timcscale at the cloud scale, the toroidal component de- 
velops very little in such a region. Although the toroidal 
field is generated around the protostar by disk rotation, 
the polarization degree is as low as <;5% at the cloud 
scale. Thus, it is difficult to determine the toroidal fields 
from such observations. 

Another way to detect the toroidal field is to investi- 
gate the direction of the polarization vectors on the po- 
larization map with a viewing angle of 9 > 0°. That is, as 
shown in the bottom panels of FigurcllO| the polarization 
vectors just above and below the protostar (r < 500 AU) 
do not point to the cloud center; they are almost perpen- 
dicular to the global field lines andparallel to the £-axis. 
In addition, the polarization degree in this region is ex- 
tremely low. This is because the foreground and back- 
ground dust alignments cancel each other. Although it is 
difficult to determine the polarization direction in such a 
low-polarization region, the existence of the region of ex- 
tremely low polarization above and below the protostar 
is an indicator of the toroidal field. 

The toroidal field can also be confirmed by observing a 
point symmetry of the low-polarization region around th e 
protostar. As described in §4.2.31 and iTomisakal (|2011[ ). 
both the viewing angle effect and the toroidal component 
of the magnetic field produce a point-symmetric distri- 
bution of a low polarization degree with respect to the 
protostar on the polarization map (compare Figures [6] 
and I10j). The point-symmetric polarization distribution 
implies the existence of a toroidal field. Although we can- 
not simply determine the existence of the toroidal field, 
we can confirm it by combining evidence appearing in 
the polarization map, such as the direction of the polar- 
ization vectors and the distribution and symmetry of the 
polarization degree. 

5.3. Comparison with Observations 

In this subsection, we compare our results with obser- 
vations. We can create a polarization map from simu- 
lation data at an arbitrary viewing angle to reproduce 
polarization observation s. First, we c o mpar e our results 
with t he observatio n s of | Girart et al.l ()2009f ). Using the 
SMA, IGirart et al.1 ((20091) " observed a molecular cloud 
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Fig. 16. — Distributions of polarization degree times column 
density (colors), which is proportional to polarized intensity; po- 
larization vector (lines); and column density (contours, in cm -2 ) 
for model 1 at a viewing angle of (8, <j>) = (60°, 0°). The proto- 
stellar mass is M ps = 0.5 Mq at this epoch. 



core in star-forming region G31.41+0.31 and showed the 
polarization dist ribution. To fit the polarization obser- 
vation shown in IGirart et al.l ([2009D . we chose model 1 
that shows a clear hourglass structure in three dimen- 
sions. Figure [TBI shows the polarization degree times the 
column density (colors) , which is proportional to the po- 
larized intensity; the polarization vector (lines); and the 
column density (contours) for a viewing a ngle of (9, <f>) 
= (60 °, 0°). The left panel of Figure 1 of IGirart et "all 
(|2009f ) shows two peaks of polarized intensity above and 
below the protostar. In addition, a low-polarization re- 
gion appears next to the protostar in the direction per- 
pendicular to the expected global magnetic lines. These 
features are reproduced well in Figure [16J Note that 
these features do not clearly appear at a viewing angle 
of 9 — 0°, as shown in Figure [6] As described in §4.2 11 
these features (two peaks and low-polarization regions) 
are emphasized as the viewing angle 9 increases. Thus, it 
is believed that this star-forming core was observed with 
9 with respect to the global field lines. Therefore, 
this structure is very sensitive to the inclination angle. 
In addition, we can imagine not a point symmetry but 
a line symmetry o f polarization intensity in Figure 1 of 
IGirart et al.l (|2009f ). This indicates that a magnetic field, 
not rotation, controls the cloud evolution because a point 
symmetry appears when the rotation energy dominates 
the magnetic energy (sec Figure I12j) . A line symmetry is 
also confirmed in Figure 1161 

Next, we fo c us on the polarization observation of 
iShinnaga et al.l (|2012t ). They measured polarized dust 
emission toward the high-mass star-forming clump IRAS 
20126+41 04. In contrast to the star-forming cloud ob- 
served by IGirart et al.l <|2009f ) . the cloud rotation seems 
to significantly affect the dynamical cloud evolution. 
IShinnaga et al.l <|2012f ) showed the iS-shaped distribution 
of the polarization vectors, which indicates that (i) the 
magnetic field lines are not aligned with the rotation 
axis and (ii) the rotational energy is larger than or com- 
parable to the magnetic energy in the cloud. Fortu- 



nately, the star-forming clump IRAS 20126+4104 was 
observed at various wavelengths, and the directions of 
the high-speed jet, mol ecular outflow, and ci rcumstcllar 
disk were determined (jShinnaga et al.ll2012f ). In addi- 
tion, the rotation axis of the entire cloud and the global 
magnetic field direction were also observed. The obser- 
vations indicate that the cloud rotation axis is inclined 
from the global magnetic field in the star-forming clump 
IRAS 20126+4104. In such a case, the magnetic field 
lines deviate from an hourglass configuration and have 
an 5-shaped configuration when the rotational energy 
dominates the magnetic energy, as described in §4.2.31 
The ff-shaped configu ration is clearly seen in Figure 1 of 
IShinnaga et al.l (|2012f h Thus, their polarization observa- 
tions agree well with our theoretical prediction. 

Structure of magnetic field different from hourglass 
shape is also o bserved in a low- m ass forming region IRAS 
16293 - 2422 (jRao et alJl20"09h . They observed the ob- 
ject with the SMA in the 345 GHz band and obtained 
maps with molecular lines such as H 13 CO + J = 4 — 3, 
SiO J = 8 - 7 and CO J = 3 - 2 in addition to the 
dust continuum polarization. Using both polarization 
and line emission observations, they obtained an inter- 
esting configuration between the magnetic field and the 
rotation. That is, in IRAS 16293 - 2422, the rotational 
axis seems perpendicular to the magnetic field around 
source A. This is similar to the iS-type configuration ob- 
served in Figure fTBT c) of model 3, in which the polariza- 
tion B-vector is likely to be perpendicular to the rota- 
tional axis in the very center (see also §5.5p . Therefore, 
we can explain the structure of magnetic field in this 
object partly by using S-type shape structure. However, 
since this object is binary, the binary system may greatly 
affect the polarization patterns. We postpone the polar- 
ization pattern expected for binary-forming model for a 
future separate paper. 

It should be noted that we should be careful about the 
filtering effect on the large scales of the interferometers. 
However, this study simply focuses on classifying polar- 
ization patterns with MHD simulations. Therefore, we 
do not consider any effects of interferometers. We should 
consider such effects for specific cases in future works. 

5.4. Effects of Sink Cells 

To investigate the long-term evolution of the magnetic 
field in star-forming cloud cores, we adopted sink cells, 
as described in §2.21 With sink cells, we did not resolve 
the high-density gas region around the protostar (or the 
center of the cloud). In the low-density gas region (or 
at large scale), the magnetic field is well coupled with 
neutrals and the magnetic field lines are accompanied 
by neutral gas motion. In principle, the gas motion in 
the low-density region is affected very little by the in- 
ner high-density gas region because the gas continues to 
fall onto the center of the cloud, and the gas motion is 
determined by the gravity as well as the local magnetic 
field, rotation, and pressure gradient force. Note that 
we correctly calculated the gravity inside the sink by 
adding the gas removed from inside the sink to the gravi- 
tational potential (see §2.2|) . Thus, we believe that we do 
not always resolve the region very close to the protostar 
(<C 10 AU) when we investigate the evolution of cloud- 
scale magnetic field lines (~ 1000 AU). In this study, al- 
though we adopted a sink radius of r sul k = 15 AU, we 
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could correctly investigate the magnetic field lines in the 
regi on of r ^> 15 AU. In addition , in our previous stud- 
ies (|Machida fc Matsum oto 2012), we confirmed that the 
magnetic field lines at the cloud scale changed very little 
even when we adopted a smaller sink radius. 

The protostellar outflow may affect the large-scale 
magnetic field structure. The molecular outflow is be- 
lieved to be driven by th e first core (jLarsonl H96l 
iMasunaga fc In utsuka 2000) before protostar formation 
(jTomis aka 2002) a nd by the circumstellar disk a fter pro- 
tostar formation (jMachida fc Matsumotoir2012D . In this 
study, because both the first core and the circumstel- 
lar disk are not resolved with sufficient resolution, no 
powerful outflow appears; only a weak outflow appears 
for model 2. However, both observations and numeri- 
cal simulations have shown that protostellar outflows are 
aligned with magnetic field lines. Because plasma mov- 
ing along the magnetic field lines never bends the lines, 
the outflows are considered not to significantly change 
the magnetic field configuration. Thus, we can safely ig- 
nore the effects of the protostellar outflow in studying 
the magnetic field configuration at the large scale. As a 
first step, we believe that our simple setting adopted in 
this study helps us to understand the global evolution of 
magnetic field lines in the collapsing gas cloud. 

In recent cloud-scale polarization observations, we 
coul d not clearly confirm the effect of protostella r out- 
flow ([GlrarriEaUllOOl HH; [ShinHagaeEaDIloil) . The 
reason is that the spatial resolution of such observations 
is not sufficient to resolve the effects of protostellar out- 
flows. When the protostellar outflow is driven by mag- 
netocentrifugal or magnetic pressure gradient forces, the 
toroidal component of the magnetic field dominates the 
p oloidal field inside the outflow. In such shown 
in iTomisakal (|2011| ) , a very low polarization degree is re- 
alized inside the outflow because the toroidal components 
canceled out each other. Thus, we must consider the ef- 
fects of protostellar outflow in order to investigate the 
polarization in the star-forming core at a small scale. To 
precisely investigate the effect of the magnetic configu- 
ration on the star formation process, we need a higher 
spatial resolution in both simulations and observations 
to resolve the protostellar outflow in future. 

5.5. Velocity Map 

The gas motion is related to the configuration of mag- 
netic field lines. For example, the magnetic field lines 
(or the polarization vectors) are highly distorted by the 
cloud rotation as seen in Figures QT] and [T2J Thus, the 
velocity information is useful to understand the config- 
uration and evolution of the magnetic field lines. To 
investigate the relation between magnetic field lines and 
velocity fields, we calculated the first moment of velocity 
along the line-of-sight, 



J pv • nds 

J pais 



(20) 



where v is velocity vector at each point and n is the unit 
vector along the line-of-sight. With the average veloc- 
ity [cq. (|20[) ]. we made the velocity map for each model. 
Figure [T7] shows the velocity distribution in the range 
of < 0.15 kms -1 for models 2, 3 and 4, in which 
each panel corresponds to the polarization map of Fig- 



ures HH and [T3] Note that the column density contour 
is superimposed for reference in Figure 1171 

Figure |T7la) is the edge-on view [((9, 0) = (90°, 0°)] 
of velocity map for model 2, and corresponds to the po- 
larization map of Figure ^c). In model 2, the initial 
rotation axis is parallel to the initial magnetic field lines. 
From this panel, we can confirm a symmetric structure 
between blueshifted and redshifted velocity contours, in- 
dicating that the rotation axis of infalling envelope and 
disk is parallel to the global magnetic field lines (or r/- 
and z-axis) . Figures I17f b) and (c) are the velocity maps 
for model 3, and correspond to the polarization maps of 
Figures I13f a) and (c), in which the viewing angles are 
(9, 0) = (90°, 0°) for Figure [TTTb) and (90°, 90°) for 
Figure IT7T c'). As described in §4.2.3, for model 3, the 
initial rotation axis is inclined from the initial magnetic 
field lines, and thus the disk normal is not parallel to 
the magnetic field lines [see Fig. fTBT b)]. In Figures fTW b) 
and (c), the black contour indicating (v) = near the 
center of the cloud (r <; 1000 AU) roughly corresponds 
to the projected disk rotation axis because the disk rota- 
tion velocity dominates the infall velocity in this region. 
In Figure IT7T b). the disk is tilted as the near side of 
the disk is lower and the far side is upper with respect 
to the line-of-sight. On the other hand, we see the disk 
from the edge in Figure IT7T c'). Thus, the rotation pattern 
in Figure \T7\c) is more clear than that in Figure IT7T b) . 
Comparison of Figure rTBT b) with Figures fT7T b) and (c) 
indicates that the disk rotation axis is parallel to the 
disk normal. Thus, we can identify the disk normal di- 
rection also from the velocity information. As described 
in §4.2.3, the magnetic field lines deviate from the hour- 
glass configuration when the disk normal is not parallel 
to the magnetic field lines. 

In contrast to models 2 and 3, the infall motion is 
more emphasized on the velocity map for model 4. Fig- 
ures [P7Td) and (e) are the velocity maps for model 4, and 
correspond to the polarization maps of Figurcs ll4f a) and 
(b), in which the viewing angles are (6, <fi) = (90°, 0°) 
for Figure E^d) and (90°, 90°) for Figure QTfe). For 
model 4, the initial rotation axis is inclined from the ini- 
tial magnetic field lines as in model 3, while the initial 
strength of magnetic field is stronger than that for model 

3. In Figures flTTd) and (e), the column density distri- 
bution implies that the disk normal is roughly parallel 
to 77-axis (or z-axis). On the other hand, the velocity 
contour of (v) = is almost perpendicular to the 77-axis 
around the center of the cloud (r < 1000 AU) and the ve- 
locity gradient is observed in the 77-direction. For model 

4, a strong magnetic field effectively transfers the angu- 
lar momentum by the magnetic braking mechanism as 
described in §.4.1. As a result, a tiny rotating disk is 
formed in the proximity of the protostar (r <C 1000 AU). 
Note that the disk-like structure of the column density 
contour in this panel corresponds to a pseudo-disk that is 
not supported by the rotation. In addition, the rotation 
of the infalling envelope is also transferred by the mag- 
netic braking. Therefore, in these panels, the rotation 
motion is not apparent; the infall motion dominates the 
rotation motion especially in the scale of ~ 1000 AU. In 
Figure fT7T d) , the disk is tilted as the near side is down 
and the far side is up. Therefore, infalling gas along the 
disk is observed as a redshifted one in r/ < region but 
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Fig. 17. — Velocity (red, blue and black contours) and column density (gray scale) maps for models 2, 3, and 4. In each panel, 
rcdshiftcd velocities of (v) = —0.15 kms -1 , — O.lkrns" 1 , — 0.075 kms -1 , — 0.05kms — 1 and —O.Olkms -1 are plotted with red contours 
while blucshiftcd velocities of (v) = O.Olkms -1 , 0.05 kms -1 , 0.075 kms -1 , O.lkrns -1 and 0.15kms — 1 are plotted with blu e co nt ours . The 
black c onto ur corresponds to (v) = 0. Panels (a), (b), (c), (d) and (e) correspond to the polarization maps of Figures l9Tc~). IT3Ta^. I13f c^. 1141 
(a) and !14f b'). respectively. 



as a blueshifted gas in rj > region. In contrast, the disk 
is tilted as the near side is up and the far side is down 
in Figure [TTJe). Thus, the infall motion is observed in 
an opposite way to (d). Therefore, with the velocity in- 
formation, we can understand the relation between the 
configuration of magnetic field and the disk or infalling 
envelope. 

6. SUMMARY 

In this study, after we calculated the evolution of 
clouds with different initial parameters including mag- 
netic field strength Bq, rotation rate Qq, and angle Sq 
between the magnetic field and the rotation axis, we cre- 
ated polarization maps from the data taken from three- 
dimensional MHD simulations. Then, we investigated 
the polarization distribution and polarization vectors at 
different epochs with different viewing angles for each 
model and obtained the following results. 

• In the prestellar stage, the polarization distribu- 
tions of the models show no apparent features. All 
the models exhibited a gradual decrease in the po- 
larization degree around the center of the collapsing 
cloud. Strictly speaking, slight differences among 
the models appeared in the polarization distribu- 
tion, depending on the initial cloud parameters of 
flo and 5q and the viewing angle. However, we ex- 
pect that such slight differences cannot be distin- 
guished even in future polarization observations. 

• After protostar formation (i.e., in the protostellar 
phase), the polarization distribution showed some 
clear differences among the models. In addition, 
it depended strongly on the viewing angle. De- 
pending on the viewing angle, a region of very low 
polarization appears next to the protostar because 



the radial component of the magnetic field around 
the protostar was canceled out. Using the low- 
polarization region next to the protostar, we can es- 
timate the angle between the direction of the global 
magnetic field and the line of sight. 

• Another difference is caused by the toroidal field. 
In the protostellar phase, cloud rotation generates 
the toroidal field, which creates a region of ex- 
tremely low polarization above and below the pro- 
tostar. On the other hand, no strong toroidal field 
appears in the prestellar stage because the freefall 
timescale is shorter than the rotation timescale. 
Thus, we can predict the existence of the embedded 
protostar and circumstellar disk from such a region 
of extremely low polarization around the protostar. 

• The other difference is the symmetry of the po- 
larization distribution. Both the rotation and the 
toroidal field break a line symmetry and produce a 
point symmetry with respect to the position of the 
protostar. The point-symmetric distribution of the 
polarization is an indicator of non-negligible rota- 
tion and a toroidal field. 

• The configuration of the polarization vectors de- 
pends strongly on the viewing angle in both the 
prestellar and protostellar stages. We found that 
an hourglass-shaped magnetic field line structure 
is not always realized in the collapsing cloud. In- 
stead, an .S'-shapcd configuration of polarization 
vectors (or magnetic field lines) often appears, de- 
pending on the viewing angle, when the global mag- 
netic field lines are not aligned with the cloud ro- 
tation axis. The S'-shapcd configuration is gradu- 
ally transformed into the hourglass configuration. 
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Thus, we observed the S'-shapcd configuration of 
the magnetic field in the early evolutionary stage. 

• The S'-shaped configuration appears in clouds hav- 
ing a rotational energy larger than the magnetic 
energy. On the other hand, the hourglass configu- 
ration always appears independent of the viewing 
angle in clouds having a magnetic energy larger 
than the rotational energy. Thus, from the config- 
uration of the polarization vectors, we can predict 
the properties of the host cloud and the evolution- 
ary stage of the protostar. 

This study complements future polarization observa- 
tions. Without this type of effort, we cannot understand 
the real evolution and configuration of the magnetic field 
because the observed magnetic field is projected onto the 
celestial plane. In this study, we showed the polariza- 
tion distribution and polarization vectors while resolv- 
ing considerably small-scale structures (~ 10 AU). Our 
results agree well with recent observations at the large 
(or cloud) scale, and we can determine their evolution- 
ary stage and cloud properties. We also created veloc- 
ity maps to compare our results with observations, and 



showed that the combination of velocity and polarization 
information gives us a better understanding of the con- 
figuration of magnetic field lines and properties of the 
collapsing cloud core. Although recent observations still 
lack sufficient spatial resolution to investigate small-scale 
structures around the protostar, our study is useful for 
understanding the effect of the magnetic field on the star 
formation process at a deeper level with near-future in- 
struments such as ALMA. 
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